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ABSTRACT 

Based on SINFONI Ha, [Nil] and [SII] AO data of 30 z - 2 star-forming galaxies (SFGs) from the 
SINS and zcSINF surveys, we find a strong correlation of the Ha broad flux fraction with the star 
formation surface density of the galaxy, with an apparent threshold for strong outflows occurring at 
1 MQyr~^ kpc~^. Above this threshold, we find that SFGs with logm,>10 have similar or perhaps 
greater wind mass loading factors (77 = Mout/SFR) and faster outflow velocities than lower mass SFGs. 
This trend suggests that the majority of outflowing gas at z ~ 2 may derive from high-mass SFGs, and 
that the z 2 mass-metallicity relation is driven more by dilution of enriched gas in the galaxy gas 
reservoir than by the efhciency of outflows. The mass loading factor is also correlated with the SFR 
and inclination, such that more star-forming and face-on galaxies launch more powerful outflows. For 
galaxies that have evidence for strong outflows, we find that the broad emission is spatially extended 
to at least the half-light radius (~ a few kpc). We propose that the observed threshold for strong 
outflows and the observed mass loading of these winds can be explained by a simple model wherein 
break-out of winds is governed by pressure balance in the disk. Using the ratio of the [SII] doublet in a 
broad and narrow component, we find that outflowing gas has a density of ^10-100 cm^'^, significantly 
less than that of the star forming gas (600 cm~'^). 

Subject headings: cosmology: observations — galaxies: high-redshift — galaxies: evolution — infrared: 
galaxies 



1. INTRODUCTION 



Based on observations at the Very Large Telescope (VLT) 
of the European Southern Observatory (ESO), Paranal, Chile 
(ESO program IDs 076.A-0527, 079.A-0341, 080.A-0330, 080.A- 
0339, 080.A-0635, 183.A-0781). 
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Most high-z SFGs show evidence for powerful galactic 
outflows, as indicat e d by UV absorption spectroscopy 
(iPettini et al.l [2OQOI: iShaplev et al.l [200l iSteidel et al.l 
120101: IWeiner et al.l 120091: IKornei et al.l 120121) and 
broad Ha emission-line profiles fShapiro. "eTaLl [20091 : 
Genzel et al. 2011; Newman et al. 2012). This 'star for- 
mation feedback' may be an essential ingredient in the 
evolution of high-z star forming galaxies. However, little 
is yet known about what drives these outflows. 

In this letter we analyze the outflow proper- 
tie s of 30 z ~ 2 SFGs , discu ssed in more detail 
in iForster-Schreiber et al.l (|2012 h and based on new 
high-quality Ha emission line SINFONI/ VLT integral 
field (IFU) spect rosco py with adaptive optics (AO) 
(|Eisenhauer et al.ll2003l:lBonnet et al.ll2004D . We adopt a 
ACDM cos mology with D.^, =0-2 7, »h = 0.046 and Hp =70 
km/s/Mpc ("Komat su et al.ll2011| ). and a lChabrieri (|2003l ) 
initial stellar mass function (IMF). 

2. OBSERVATIONS, DATA REDUCTION AND 
ANALYSIS TECHNIQUES 

The galaxy sampl e discussed here is desc ri bed i n 
iMancini et all (|2011|) : iForster-Schreiber et al.l (|2012D . 
Our 30 z~2-2.5 SFGs are drawn from the SINS and 
zC-SINF surveys of Ha-f [Nil] integral field spectr oscopy 
with SINFONI on th e ESO VLT (Forstcr Schreiber et alj 
l2009HMancim et al.llM l: Forste r-Schreiber et al. 2012^ 
They sample the z~2.2 star-forming 'main sequence' for 
stellar masses between 10^'^ to lO"'^"'^'^ Mq. Our sample 
includes 27 galaxies observed with AO-mode with suf- 
ficient SNR for analysis (0.05" pixels, typical FWHM 
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0.2-0.25") in addition to 3 galaxies observed in seeing 
limited mode that have very extended disks and thus 
are well resolved with 0.125" pixels and FWHM~0.5". 
The data are of high quality (~ 6h average integration 
time per galaxy, with a range from 2 to 22h, and a to- 
tal integration time for the entire sample of 180h) and 
were reduced with o ur standard data, reduction methods 
and analysis tools (iSchreiber et al.l 120041: iPaviesI [20071: 
iForster Schreiber et al.ll2009l: iMancini et a l. 2011). 

Stellar masses and ages are derived from SED modeling 
inHlrster Schr eiber et al. (200^; IMancini etall ()2011l ) 
assuming constan t or exponentially declinin g star forma- 
tion histories with lBruzual fc Charlotl ()2003l ) tracks. The 
SFRs are derived from the SED mode li ng an d from Ha 
(SF R=Lff»/2 . 1 x10'^^ erg/s, iKennicuttl (I1998D. co r rected 
for a lChabriei< (l200l IMF) with a lCalzetti et al.l (l2000h 
reddening law with Ay.gas = 2.3 x Ay^sED- Half light 
radii (R1/2) were computed from the cumulative flux pro- 
file determined after fitting a 2D exponential (Sersic n — 
1) to the data. We calculate molecular gas masses (Mgas) 
and surface densities {T,gas— 0.5 x fron^ 
the Ha-derived SFRs using the molecular gas to star 
formation surface density relation with ^I^s(Mq) = 
5.8x10^ X SFR (Moyr-i) (iTacconi et al.l[20B ). 

To isolate the broad emission for coadding, we re- 
move large-scale rotational velocity gradients from the 
data cubes and create spatially-integrated spectr a (see: 
IForster Schreiber et al.ll2009HShapiro et al.ll2009( ). Vari- 
ous stacks of spectra are generated by removing a fitted 
continuum from each spectrum, placing the spectra on 
a common rest-frame wavelength axis, weighting each 
spectrum by the signal to noise ratio (SNR) of Ha , and 
coadding the final spectra. 

Broad emission has previo usly been det ec ted i n 
SINS/zC-SIN F gal ax ies by 'Shapiro et alj ([20M ): 
iGenzel et eUl (|2011[ ): iNewman et al. (2012) and at- 
tributed to galactic outflows. We quantify the fraction 
of the emission line flux in this underlying broad compo- 
nent by simultaneously fitting two Gaussian functions to 
each of the Ha, [NII]AA6548,6584 and [SII]AA6716,6731 
emission lines, setting the kinematics (velocity, line 
width) of all of the narrow Gaussian components to be 
equal, and likewise with the broad components. The 
broad flux fraction is the ratio of the total flux in the 
broad emission line components (one for each emission 
line) to the total flux in all of the narrow emission line 
components. 

We also generate position-velocity (pv) data by ex- 
tracting the spectra along the major galaxy axis for each 
shifted cube with a slit wide enough to cover most of the 
minor axis emission. The data are then co-added and 
re-sampled onto a normalized spatial scale of 0.07xRi/2 
(FWHM - 0.3xRi/2) and spectral scale of 40 km/s/pixel 
(-90 km/s FWHM). 
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We ascertain that the outflows are driven by the star 
formation activity and not by an AGN as we have ex- 
cluded any known AGN from our sample, find that the 
broad emission is extended and not concentrated towards 
the kinematic or morphological center (see section 3.1), 
and measure [Nil] /Ha ratios compatible with the ex- 
pected metallicity based on the z ~ 2 mass- met allicity 




2000 3000 



Fig. 1. — Co-added (pv) diagram for the 10 sources with the best 
evidence for broad emission, with the vertical (spatial) axis nor- 
malized to Rx/2 (lower panel) and spectra from the inner (R<0.5 x 
Rl/2i blue, solid line) and outer (0.5 x Ri/2<R<R.i/2i red, dashed 
lino) regions of the galaxies (upper panel). The black contours 
show the unnormalized stack, while the colors have been normal- 
ized according to the peak Ha flux for each spatial row. The broad 
emission in light blue (> 200-300 km/s from the systemic velocity) 
is spatially extended (vertical axis) out to at least Ri/2- The inner 
and outer region spectra have Fj^o^^/F^nrrom = 1-21 it 0.039 and 
0.89 ± 0.049, respectively. 

relation ()Erb et al.l I2006D and too low to be consistent 
with the presence of an AGN. In several cases we trace 
the origin of the outflows to individual g iant star form- 
ing clumps enibedd ed in the disk (Gcnzel et al.l 120111 : 
INewman et al.ll2012f ). 

3.1. Outflows are Spatially Extended 

The broad emission (FWHM '--^500 km/s) is spatially 
extended over the half-light radii (R1/2), which corre- 
sponds to a few kpc in the small galaxies and up to 6 
kpc in the larger disks. This is demonstrated in the bot- 
tom panel of Figure 1, where we show a co- added pv 
diagram of the line emission (large-scale velocity gradi- 
ents removed) along the morphological major axis for 
the galaxies. Broad emission is detected at least out 
to Ri/2, and its width between 0.5 x R1/2 and R1/2 
(FWHM = 550 km/s) is nearly as broad as within 0.5 
X Ri/2 (FWHM = 610 km/s). The upper panel of Fig- 
ure 1 shows that the emission in the outer regions of 
the galaxies is almost as broad as in the inner regions, 
and implies a constant or decreasing mass-loading with 
increasing galacto-centric radius. 

3.2. The Broad Flux Fraction is Strongly Dependent 
on the SF Surface Density 

We explore the dependence of the broad flux fraction 
on galaxy properties by dividing the entire sample into 
two bins each (low/high) in star formation rate, stellar 
mass, size, inclination and star formation surface density. 
We then co-added the spectra in each of the bins and 
computed the ratio of the broad to narrow Ha -|-nebular 
emission. The results are summarized in Table 1. 

The star formation surface density has the largest 
effect on the broad flux fraction, such that galaxies 
with Y}sFR>^ M0j/r~^fcpc~^ drive the strongest out- 
flows. This result is consistent with what was observed 
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TABLE 1 

Broad flux fraction for high and low bins of different galaxy properties 

Property Abroad/ Fnarrow FWHMbroad (km/s) A Vbroad - narrow (km/s) Significance Dividing 

High bin Low bin High bin Low bin High bin Low bin of difference value 



SFR 


0.65 ± 0.074 


0.47 ± 0.042 


510 


423 


-23 


-23 


2a 


100 Me/yr 


R-l/2 


0.50 ± 0.054 


0.76 ± 0.082 


503 


435 


-15 


-35 


3(T 


3 kpc 


SsFH 


0.77 ± 0.027 


0.16 ± 0.030 


498 


423 


-25 


-43 


20<T 


1 MQyr~^fcpc~^ 


inclination 


0.47 ± 0.055 


0.77 ± 0.091 


510 


515 


-39 


-32 


3<T 


49-55° 


m. 1 


0.63 ± 0.056 


0.41 ± 0.042 


524 


423 


-29 


-25 


4cr 


1x10^° Mq 


m. = 


0.23 ± 0.21 


0.12 ± 0.040 


423 


423 


-28 


-48 


0.5cr 


m. 3 


0.85 ± 0.097 


0.68 ± 0.10 


522 


435 


-25 


-14 


1.5cr 





The high/low bins for SFR, m*, ^sfr, Ri/25 a^nd inclination are divided above and below the value(s) shown in column 9 and have roughly an 
equal number of galaxies in each bin. The spectrum from each stack was fit constraining the narrow FWHM to between 190 and 210 km/s, so that 
the broad component was not fit as a very broad narrow component, and allowing the relative velocities of the two components and the broad line 
width (FWHM) to vary. The best fit broad FWHM and the relative velocities for each bin arc shown in columns 4—7. For the m* bins, (1) is for 
all galaxies, (2) is for galaxies with T1sfr<^ yiQyr'^^kpc'^'^ ^ and (3) is for galaxies with SsFf£>l MQyr ~ "'" fcpc~^ . 



in z - SFGs bv lChen et all ()2010D . where they found a 
strong correlation of Na D absorption equivalent width 
(from SDSS data of ~ 100,000 galaxies) with T^sfr- 

In Figure 2, we see the dependence of the broad Ha 
flux fraction on Ss^^from 5 S5^fl-binned points as 
well as three massive star-forming clumps. The trend 
seen in Figure 2 can be well described as a 'thresh- 
old' for outflows. All of the data points with Ssi="fl >1 
MQyr~^kpc~^ have F^road/F narrow > 0.7, while all those 
below this threshold have F broad /F narrow < 0.25. 

The broad flux fraction also varies (to a lesser extent) 
with m^,, SFR, R1/2 and inclination such that stronger 
outflows are observed for massive, high SFR, compact, 
face-on galaxies. We further explore the m* trend by di- 
viding each of the high and low EsFfi bins into high and 
low m* bins. Unsurprisingly, we find that both of the m* 
bins below this threshold show negligible evidence for 
broad emission, and both high I]5i?flbins show strong 
outflows. Above the T,sfr threshold, the broad flux frac- 
tion is similar or perhaps somewhat larger in the high m, 
bin. The FWHM of the broad line is larger for the high 
m, bin (522 km/s vs. 435 km/s) and more blueshifted 
(-25 km/s vs. -14 km/s). This increase is smaller than a 
'viria l' scaling, Vn,, * < Vesc ~ Vc ~ m*°-^'^ ()Murrav et al.l 
120051: iMartinI 120051 ). which would imply a factor of ~ 1.8 
in Vout between the two mass bins (logm,= 9.85 and 
10.61). 

The SFR trend is linked to that of m*, due to the 
strong correlation of m, and SFR for 'main-sequence' 
SFGs at all redshifts. The R1/2 trend is likely due to the 
higher Sgi^jj values for the more compact galaxies. The 
trend with inclination supports a bipolar outflow emerg- 
ing perpendicularly from the pla ne of the star-forrn ing 
disk, c onsistent wit h the results of lKor nei et al.' (2012*, at 
D.lChen et all (2010. at z ~ O UBouche et ah (2012 1 



J2012 



and lBordoloi et alj |201lf ) (but see lLaw et al.l 1)201 
3.3. Local Electron Density of the Outflow 

We estimate the ratio of the [SII] doublet in the broad 
and narrow components for a stack of the galaxy spec- 
tra, in order to constrain the star-forming gas and wind 
densities. The 14 galaxies used in this stack are selected 
such that they do not have strong OH sky features close 
to the location of the [SII] lines and have noticeable broad 
Ha components. 

For the narrow component, we find 
F([SII]A6716)/F([SII]A6731) = 0.99 ± 0.27, and for 
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Fig. 2. — Dependence of the broad to narrow Ha flux ratio on 
12$ PH. Above: Stacks in 5 star formation surface density bins with 
5-8 galaxies each. The right vertical axis shows the corresponding 
mass loading with the model explained in the text and wind param- 
eters listed in the figure. The black points are the stacked galaxies 
and the colored symbols are individual star-forming clumps. The 
horizontal error bars for the binned points represent the range in 
^SFR for the galaxies in each bin. The vertical error bars represent 
the range in Fi„.aad/Pnarrow values achieved by varying the stacks 
(for the binned points) and the measurement error (for the clumps). 
Below; Co-added Ha -(-[Nil] spectra of low and high T,sfr bins 
with the data in black, the best fit shown in red, the narrow com- 
ponents shown in blue, and the broad components shown in green 
with grey shading. 



the broad component, we find that the ratio = 1.43 ± 
0.40, corresponding to electron densities of 600 (-1-1000/- 
450) and 10 (-f590/-10) cm~^, respectively (|Osterbrockl 
Il989f ). Although our estimate is uncertain, the inferred 
local gas density in the outflow is plausibly more than 
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0.650 0.655 0.660 0.655 0.670 0.675 

Wavelength (^m) 

Fig. 3. — Co-added spectrum of 14 SFGs with [SII] detections 
that are free of OH contamination. The black Une is the data, 
the red line is the best fit, and the green lines with grey shading 
(with FWHM = 407 km/s) and blue dashed lines (with FWHM = 
197 km/s) (seen in the inset) are the broad and narrow Gaussian 
fits, respectively. We follow the fitting method described in the 
text, except here we allow the amplitudes of the broad and narrow 
components in all nebular lines to vary. 

an order of magnitude less than the density of the 
star- forming regions (traced by the narrow component), 
consistent with a diffuse outflow. Assuming th e wind 
density scaling law presented in ICecil et all (|2001[ ) 
from observations of M82 (ug — 5.7 x where f 

is the dense gas filling factor of the wind), our wind 
density implies a filling factor of 50%. As this filling 
factor is perhaps too high, our wind density is likely a 
lower limit , ownin g to the low SNR of the [SII] lines. 
ICecil et al.l (I2OOII) also find a minimum filling factor 
of 0.3%, implying a wind density of 100 cm~^. Thus 
we take an average of these lower and upper limits, 
resulting in a wind density of 50 cm~^. Figure 3 shows 
the co-added spectrum along with the best fit. 

3.4. Mass-loading of High-z Galactic Outflows 

We assum e the simp lified outflow mode l of 
IGenzel et all (poTll ) and INewman etHI (poTl ) for 

a warm ionized outflow with a radially constant 
outflow velocity and mass loss rate, to convert the 
FbroadI Fnarrow ratio in Figure 2 into a mass-loading 
factor ry = Mout/SFR, 



Mo 



Mw X 



Vput 
Rnn.t 



l.36mH 

iHane 



Lhcx X 



Fbr. 



oad 



Vout 



Ftotal ) Rout 
(1) 

where Mout is the mass outflow rate, M^t, is the instan- 
taneous mass in the outflow, Vout is the average veloc- 
ity of the outflow. Rout is the radial extent, tclh is the 
atomic mass, ^Ha is the Ha emissivity at Te = IC* K 
{iHa = 3.56 X 10~^^ erg cm~'^ s~^), Ue is the local elec- 
tron density in the outflow, Lna is the total extinction- 
corrected Ha luminosity, and Fbroad/Ftotai is the fraction 
of the total fiux in the broad component. We assume Vout 
= 400 km/s (see section 4) and Rout = 3 kpc. The latter 
is based on the spatial offset of outflow ing gas from a 
massiv e star-forming clump described in INewman et al.l 
(|2012f) . and from the flnding in section 3.1 that the broad 
emission is extended to several kpc from the galaxy cen- 
ter. The emission line based estimate is independent of 
the coUimation of the outflow. We adopt Ug = 50 cm~'^ 



as the average local electron density, derived in the pre- 
vious section. 

The inferred mass loading factors corresponding to 
the broad flux fraction are shown on the right-most y- 
axis in Figure 2. The mass- loading factors for galax- 
ies below the E^i^fl threshold are ^ 0.5 and above the 
threshold are ~ 2, albeit with an overall absolute un- 
certainty of at least a factor of 3. Despite the uncer- 
tainties, a mass loading factor of 2 is consistent with ob- 
servations of both local starbursts ("Heck man et al.lFlQQOl : 
lyeflleux et al 2005; Chen et al. 2010) and high-z SFGs 



120101 IGenzel et al.l 20111 Bouche et al. 


2612I), as well as 


theoretical predictions (iMurrav et al.l 


20051 iDave et al.l 


l201lHHoDkins et al.ll2012D. 



4. WHY DO THE OUTFLOWS DEPEND SO 
STRONGLY ON T^sfr'^ 

We propose that the strong dependence of 77 on T,sfr 
discussed in section 3.4 is mainly caused by the threshold 
that governs when star formation feedback can break out 
of the dense gas layer in t he disk. 

Following IMurrav et all ([2005. ) and lOstriker fc Shettvl 
in a baryon dominated galactic disk in pressure 
equilibrium, the weight of the disk balances the pressure 
generated from star-formation feedback in the form of 
supernovae, stellar winds, HII regions, cosmic rays and 
radiation from the star-forming layer. If that pressure 
exceeds the weight of the disk, then a momentum-driven 
outflow is launched perpendicular to the galactic plane, 
with a threshold T,sfr {MQyr^^ kpc~^) of. 



^SFR,th 



^Gfg ^2 = 9 X ■^3-°-^ ^ ^^.500 

2{Ptot/m^) (ftot/m,) 1000 



(2) 



where G is the gravitational constant, (Ptot/m*) 
is the characteristic total momentum injection per 
mass, Ed is the disk surface density, and fg is 
the gas fraction, with fiducial values of 1000 km/s, 
500 MfT^pc"^ and 0.5, respectively (Erb et al.' '2007 
Tacconi et alEo lO: Daddi et al. 2010; Ostriker & Shett 
20111: iForster Schreiber et al.l 1200911 . This is simflar 
to the Eddin g ton h mit of momentum driven winds 
(jMurrav et al.l |2005[) . which gives ^sFR,threshoid oc 
v'^R~^, where Vc is the circular velocity. 
The wind has a total momentum outflow rate. 



Mw-foo < 



p, 



tot 



SFR 



(3) 



where Vqo is the outflow velocity at large distances from 
the mid-plane, implying a mass-loading factor, 



77 



SFR 



< 



(Ptot/m^) 



2.5 X 



(-Ptot /to*) 1000 
foo.400 



(4) 



with Voo,4oo in units of 400 km/s. This velocity is moti- 
vated by the observed line widths of broad emission-line 
components in both l ow-z and high-z SFGs (Pettini et aQ 
2000 1 : iMartinI [2001 iVeflleux etall [20051: iWeiner et all 



2003: iSteidel et al.M2010l: IGenzel et al.l 120111) as well as 
from this work. The adopted value of the characteristic 
momentum injection is based on work by iMurrav et ahl 
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((20051) : IVeilleux et aP (12001 : lOstriker fc Shettvl ([20lTh 
wherein the direct momentum injection contributions of 
supernovae, steflar winds, HII regions and radiation pres- 
sure are afl roughly comparable, with (P/m») ^ 200-300 
km/s, giving a total (P/m*) of 1000 km/s. 

This simple model predicts that for constant 77 the 
ability of an outflow to break out of the disk strongly 
depends on the star formation surface density, and this 
threshold value scales linearly with gas fraction and as 
the square of the disk surface density. Second, the mass 
loading of galactic winds above breakout and the thresh- 
old star formation surface density for breakout both de- 
pend on (Ptot/m*). Thus, if one of these quantities can 
be determined for a galaxy of known properties, then 
the other can be predicted. The model also predicts 
that it is harder to launch a wind from a more mas- 
sive galaxy (with larger v^/R^). This final prediction is 
also acquired by writing the threshold i n terms of m,, 
by assuming (1) the Ri/2-m* relation of llchikawa et al.l 
((Ml (R1/2-- m,o i4), (2) the SFR-m* relation for nor- 
mal SFGs (SFR = m*0-^ x ((l-Hz)/3. 2)2-6) and (3) a con- 
stant gas depletion timescale {idepi = Mg/SFR 5.8 x 
10^ vr. iTacconi etall (HoH)). This gives an m* thresh- 
old of few X 10"^° Mq, such that it is easier for outflows 
to 'break-out' from galaxies below this threshold. 

The first two predictions are met by the avail- 
able data. For z^O normal SFGs with fg ~0.07 
and Ed ^500 Mopc"^, Equation 4 suggests a criti- 
cal break-out star formation surface density of ~0.1 
M0yr~^fcpc~2. And indeed, winds are only detected 
in compact star forming dwarfs and/or starbursts above 
~0.1 M(7)j!/r~^fciPc-2(lHeckman et al.lll990HVeilleux et al.l 
[2OO5I: IChen et al.l [20101 ). see also IMurrav et al.l (|201in . 
Mass loading factors for galaxies above this threshold 
are estimated to be around 1. 

For the typical high-z SFGs presented in this paper 
(fg - 0.5, - 500-1000 Mopc-2)^ Equation 2 sug- 
gests a breakout star formation surface density near or 
slightly above 1 MQi/r^^fcpc"^, in remarkable agreement 
with our findings in section 3.2 and Figure 2. Above this 
threshold. Equation 4 predicts a mass loading of ^2.5, 
again in good agreement with our observations. 

A positive co rrelation of 77 wi th Yisfr (which was 
also observed bv lChen et al.l (|2010t )) is possible evidence 
against the scenario in which winds are launched by the 
energy of hot supernovae bubbles, as in this case the 

— 1/2 

shorter cooling time of the dense gas suggests t] cx 'S'sfr 
((Hopkins et al.|[20T2l ). and could support the momentum- 
driven wind model. 

However, the third prediction of the simple model 
(massive galaxies are less efficient at driving winds) does 
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not appear to be met by our data. For the SFGs 
in the 'wind regime' (SsFi?, > 1 M0?/r~^fcpc~2) we do 
not see any significant variation of ^broad/^narrow as a 
function of Vc, which is expected t heoretically for mo- 
mentum driven winds (v v7 ■ \. IMurrav et al l [20051 : 
iQppenheimer fc Davell 2006l[2008t [Hopkins et al.ll2012| ). 

This finding as well as the approximate m* indepen- 
dence of the observed F broad /F narrow ratio surprisingly 
suggest that massive SFGs have mass loading factors that 
are similar to or higher than those of low mass SFGs. 
If so, the volume averaged outflow rate and metal en- 
richment of the circum-galactic/inter- galactic medium at 
z ^ 2 may be dominated by massive galaxies just around 
the Schechter mass (logM* ^ 10.65). Indeed, with a 
slope a of the main-sequence of SFGs (SFR cx m") ~ 
0.7 or higher, the SFR (and thus the outflow rate, with 
constant 77) increases faster with m* than the volume 
densit y of SFGs decreases ($(m,) ^ m^~°-^. [Peng et ahl 
((20T0l )). We note that a massive galaxy dominance in 
the metal enrichment of the intracluster medium (ICM) 
is demanded by the fact that in clusters ^ 2/3 of the 
metal mass is contained in the ICM whereas only ~ 1/3 
is still locked into stars and galaxies (Rcnzini ,1997[ ). 

Further, the z 2 m* dependence of the mass- 
metallicity relation (Z ~ m*'^'^) may be primarily driven 
by the larger (diluting) gas fractions in low m* SFGs, 
besides these galaxies driving winds more eff ectively 
(Mg/m^— tdepi X SFR/m*'^ m,~°-^), see also E rb et al.l 
(|2009D . In addition, if high-m, galaxies are indeed eject- 
ing this much mass, we may have uncovered a mechanism 
contributing to t he quenching of s tar formation near the 
Schechter mass ((Peng et al.ll2010[ ). Using simple scaling 
arguments and the median SFRs and gas masses for the 
high and low m* bins, we flnd that in an open box model 
(with the mass accretion rate = 2 x SFR and where half 
of the mass formed in stars is returned back to the ISM), 
high-m* galaxies will run out of gas in 750 Myr, whereas 
low-m* galaxies will not run out of gas for 3.3 Gyr. 

From equation 4, we would assume that if rj is con- 
stant for both m* bins, but the outflow velocity is greater 
for the high-m, galaxies, then (Ptot/m,) would have to 
be larger for these galaxies as well. Perhaps this is due 
to the higher metallicity of the high-m* galaxies, which 
would increase the dust content, potentially raising the 
radiation pressure contribution to (Ptot/ni*). 
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